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ABSTRACT 

The  upper  most  layers  of  the  ocean,  along  with  the  lower  atmospheric  boundary 
layer,  play  a  crucial  role  in  the  air-sea  fluxes  of  momentum,  heat,  and  mass,  thereby 
providing  important  boundary  conditions  for  both  the  atmosphere  and  the  oceans 
that  control  the  evolution  of  weather  and  climate.  In  particular,  the  fluxes  of  heat 
and  gas  rely  on  exchange  processes  through  the  molecular  layers,  which  are  usu¬ 
ally  located  within  the  viscous  layer,  which  is  in  turn  modulated  by  the  waves 
and  the  turbulence  at  the  free  surface.  The  understanding  of  the  multiple  interac¬ 
tions  between,  molecular  layers,  viscous  layers,  waves  and  turbulence  is  therefore 
paramount  to  an  adequate  parameterization  of  these  fluxes.  In  this  paper,  the  au¬ 
thors  present  evidence  of  a  clear  coupling  between  the  surface  waves  and  the  sur¬ 
face  turbulence.  When  averaged  over  time  scales  longer  than  the  wave  period,  this 
coupling  yields  a  spatial  relationship  between  surface  temperature,  divergence  and 
vorticity  fields  that  is  consistent  with  spatial  patterns  of  Langmuir  turbulence.  The 
resulting  surface  velocity  field  is  hyperbolic  suggesting  that  significant  stretching 
takes  place  in  the  surface  layers.  On  time  scales  for  which  the  surface  wave  field 
is  resolved,  we  show  that  the  surface  turbulence  is  modulated  by  the  waves  in  a 
manner  qualitatively  consistent  with  rapid  distortion  theory. 


February  2008 


F.  Veron,  W.K.  Melville  and  L.  Lenain 


3 


1.  Introduction 

The  upper  most  layers  of  the  ocean,  along  with  the  lower  atmospheric  boundary  layer,  play 
a  crucial  role  in  the  air-sea  fluxes  of  momentum,  heat,  and  mass,  thereby  providing  important 
boundary  conditions  for  both  the  atmosphere  and  the  oceans  that  control  the  evolution  of  weather 
and  climate. 

The  first  models  of  boundary  layers  on  both  sides  of  the  air-sea  interface  were  developed  from 
our  understanding  of  the  turbulent  flow  over  rigid  flat  surfaces  and  extended  to  the  field  after  the 
landmark  Kansas  experiment  on  the  terrestrial  boundary  layer  (Businger  et  al.,  1971).  Conse¬ 
quently,  models  of  neutrally  stratified  flows  are  based  on  the  well-known  “law  of  the  wall”  which 
depends  on  the  assumptions  of  a  constant  stress  layer  and  horizontal  homogeneity.  The  veloc¬ 
ity  profile  then  depends  on  the  stress,  the  vertical  distance  from  the  boundary,  and  the  roughness 
length,  Z(),  which  parameterizes  the  influence  of  the  roughness  elements  at  the  surface  on  the  kine¬ 
matics  and  dynamics  of  the  flow.  When  the  flow  is  stratified,  “law  of  the  wall”  similarity  theory 
is  further  modified  to  include  buoyancy  effects,  which  are  parameterized  by  the  so  called  Monin- 
Obukhov  length,  the  height  at  which  buoyancy  and  shear  production  of  turbulence  are  of  equal 
magnitude. 

The  air-sea  interface,  however,  is  dynamically  quite  different  from  a  solid  flat  surface.  For  one 
thing,  the  velocity  field  is  not  required  to  vanish  at  the  interface.  The  ocean  surface  responds  with 
drift  currents,  surface  waves  and  turbulent  eddies  over  a  broad  range  of  scales.  There  are  other 
phenomena  such  as  bubble  injection,  spray  ejection,  rainfall,  foam  and  surfactants,  which  further 
affect  the  dynamics  and  complicate  the  problem.  Consequently,  one  would  expect  the  dynamics  of 
such  an  interfacial  layer  to  be  significantly  different  from  that  over  a  solid  flat  surface  under  similar 
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forcing  conditions.  Indeed,  while  there  is  evidence  that,  in  general,  and  on  average,  the  Monin- 
Obukov  similarity  theory  holds  over  the  ocean  (Edson  and  Fairall,  1998;  Edson  et  al.,  2004), 
there  are  also  some  notable  differences.  Over  the  last  decade  or  so  it  has  also  become  apparent 
that  surface  wave  processes  can  play  an  important  role  in  the  kinematics  and  dynamics  of  the 
boundary  layers  (e.g.  Janssen,  1989,  1999;  Komen  et  al.,  1994;  Belcher  and  Hunt,  1998;  Hristov 
et  al.,  1998;  Edson  and  Fairall,  1998,  Uz  et  al.,  2002;  Sullivan  et  al.,  2004,  2007),  and  recent 
measurements  and  models  of  the  drag  of  the  sea  surface  on  the  atmosphere  at  moderate  to  high 
wind  speeds  in  fact  suggest  that  much  of  the  momentum  transfer  at  the  surface  is  supported  by  the 
form  drag  of  the  waves.  There  is  also  evidence  that  the  air-sea  heat  flux  can  be  modulated  by  the 
wavy  surface  (Sullivan  and  McWilliams,  2002;  Veron  et  al.,  2008a).  In  addition,  breaking  waves 
generate  turbulence  (Rapp  and  Melville,  1990;  Agrawal  et  al.,  1992;  Thorpe,  1993;  Melville, 
1994;  Anis  and  Mourn,  1995;  Melville,  1996;  Terray  et  al.,  1996;  Veron  and  Melville,  1999a), 
which  along  with  small  scale  Langmuir  circulations  and  coherent  structures  (Melville  et  al.,  1998; 
Veron  and  Melville,  1999b;  McWilliams  et  al.,  1997;  Veron  and  Melville,  2001;  Sullivan  et  al., 
2004,  2007)  may  lead  to  enhanced  dissipation  and  mixing  with  significant  departures  from  the 
“law  of  the  wall”  and  may  also  result  in  increased  heat  and  gas  transfer  (Jahne  et  al.,  1987;  Hasse, 
1990;  Jahne  and  HauBecker,  1998;  Zappa  et  al.,  2001;  Garbe  et  al.,  2004;  Schimpf  et  al.,  2004; 
Turney  et  al.,  2005). 

With  such  richness  in  the  phenomena  and  dynamics  at  the  surface,  one  can  also  expect  sig¬ 
nificant  interactions  between  currents,  surface  waves,  and  turbulence.  A  well  known  example 
is  Langmuir  circulation  that  results  from  the  interactions  between  vorticity  and  the  Stokes  drift 
generated  by  the  surface  waves  (Leibovich,  1983;  Thorpe  2004).  On  smaller  time  scales,  there 
is  also  evidence  that  the  turbulence  can  be  significantly  coupled  with  the  surface  waves  (Lumley 
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and  Terray,  1983;  Cheung  and  Street,  1988a,  1988b;  Thais  and  Magnaudet,  1995;  Teixeira  and 
Belcher,  2002). 

Since  the  air-sea  transfers  of  heat  and  gas  rely  in  part  on  exchanges  through  the  diffusive 
surface  layers,  and  since  these  layers  are  typically  smaller  than  the  viscous  sub-layer,  our  ability 
to  adequately  quantify  these  fluxes  depends  on  our  understanding  of  the  small-scale  turbulence 
that  controls  the  dynamics  of  the  molecular  layers1.  In  turn,  our  understanding  of  the  small  scale 
turbulence  depends  on  our  understanding  of  the  multiple  interactions  between  the  turbulence, 
currents,  and  surface  waves. 

In  this  paper,  we  present  results  from  field  experiments  that  show  evidence  of  coupling  be¬ 
tween  the  surface  temperature,  the  surface  turbulent  velocity  fields,  and  the  surface  waves.  We 
also  show  evidence  of  coupling  between  the  surface  waves  and  the  turbulent  fields.  In  section  2, 
we  summarize  the  experimental  setup.  In  section  3,  we  present  measurements  of  surface  waves 
and  the  velocity  fields  which  provide  the  basis  for  the  results  presented  here.  Section  4  shows 
the  coupling  of  the  waves  with  the  surface  kinematic  fields  and  surface  temperature  over  long 
time  scales,  when  the  wave  effects  are  averaged  in  a  manner  similar  to  that  which  leads  to  the 
Craik-Leibovich  mechanism  (CLII)  for  the  formation  of  Langmuir  circulation  (Leibovich  1983). 
Finally,  section  5  shows  evidence  of  distortion  of  the  surface  turbulence  by  the  waves  over  short 
times  scale  when  the  wave  field  is  resolved. 

2.  Experiments 

The  measurements  described  here  were  obtained  from  a  field  experiment  conducted  from  R/P 
'in  the  case  of  gas  transfer,  these  molecular  layers  can  also  be  at  the  surface  of  entrained  bubbles. 
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FLIP ,  moored  approximately  150  miles  off  the  coast  of  southern  California,  west  of  Tanner  Bank 
(32°  40. 20' IV,  119°  19.46W,  300  m  depth),  during  20-26  August  2003. 

The  main  instruments  comprised  an  integrated  active  and  passive  infrared  imaging  and  altime¬ 
try  system  (Veron  et  ah,  2008b),  and  an  eddy  covariance  atmospheric  flux  package.  Both  systems 
are  described  in  more  detail  below. 

a.  Infrared  imaging  and  altimetry  system 

The  active  and  passive  infrared  imaging  and  altimetry  system  includes  an  infrared  camera 
(Amber  Gallileo),  a  60  IV  air-cooled  CO-2  laser  (Synrad  Firestar  T60)  equipped  with  an  industrial 
marking  head  (Synrad  FH  index)  with  two  computer-controlled  galvanometers,  a  laser  altimeter 
(Riegl  LD90-3100-EHS),  a  video  camera  (Pulnix  TM-9701),  a  6  degree-of- freedom  motion  pack¬ 
age  (Watson  Gyro  E604),  and  a  single  board  computer  (PC  Pentium  4).  All  instruments  were 
enclosed  in  a  weatherproof,  air-conditioned,  aluminum  housing.  All  instruments  and  computers 
were  synchronized  together  to  within  2  ms  and  also  to  GPS  time.  The  infrared  camera  was  set  to 
record  temperature  images  (256  x  256  pixels)  at  60  Hz,  with  a  2  ms  integration  time,  yielding 
better  than  15  mK  resolution.  The  footprint  of  the  infrared  camera  was  approximately  2  m  by 
2  m.  The  video  camera  (768  x  484  pixels)  was  synchronized  to  the  infrared  camera  and  acquired 
full  frames  at  30  Hz.  The  foot  print  of  the  infrared  camera  was  contained  within  the  foot  print  of 
the  video  camera.  The  infrared  CO2  laser  and  accompanying  marking  head  were  used  to  actively 
lay  down  patterns  of  thermal  markers  on  the  ocean  surface  in  order  to  study  the  rate  of  decay  of 
an  imposed  surface  temperature  perturbation  while  tracking  the  Lagrangian  velocity,  shear,  and 
vorticity  at  the  surface.  Finally,  the  laser  altimeter  measured  the  distance  to  the  water  surface  at 
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12  kHz  (averaged  down  to  50  Hz)  with  a  foot  print  of  5  cm  diameter,  contained  within  both 
the  infrared  and  video  images.  The  system,  among  other  things,  yields  the  velocity  at  the  surface 
by  tracking  active  thermal  markers  laid  down  with  the  CO2  laser  (Thermal  Marking  Velocimetry, 
TMV)  or  by  performing  cross-correlation  analysis  on  the  passive  surface  temperature  fields  (Parti¬ 
cle  Image  Velocimetry,  PIV).  The  detailed  performance  of  the  passive  and  active  IR  measurement 
system  for  ocean-surface  kinematics  is  described  in  Veron  et  ah,  2008b. 


Fig.  1 


b.  Eddy  covariance  system 


In  addition  to  the  optical  infrared  system,  we  used  an  eddy  covariance  system  to  acquire 
supporting  meteorological  and  atmospheric  boundary  layer  flux  data.  The  eddy  covariance  sys¬ 
tem  included  a  three-axis  anemometer/thermometer  (Campbell  CSAT  3),  an  open  path  infrared 
hygro meter/ C'02  sensor  (Licor  7500),  a  relative  humidity/temperature  sensor  (Vaisala  HMP45), 
and  a  net  radiometer  (CNR1).  Turbulent  fluxes  of  momentum,  heat  and  moisture  were  calcu¬ 
lated  using  the  covariance  method  over  30-min  averages.  The  sonic  temperature  was  corrected  for 
humidity  and  pressure,  Rotation  angles  for  correcting  the  orientation  of  the  anemometer  were  ob¬ 
tained  from  the  30-min  averages  of  the  velocity  components,  and  the  latent  heat  flux  was  corrected 
for  density  variations  (Webb  et  al.,  1980).  For  the  purposes  of  this  paper,  the  good  agreement  be¬ 
tween  the  flux  covariance  data  and  estimates  using  the  Toga  Coare  3.0  algorithm  (Fairall  et  al., 
1996;  Fairall  et  al.,  2003)  supports  the  use  of  the  Monin-Obukhov  similarity  theory  for  these  open 
ocean  conditions  (Edson  et  al.,  2004)  and  gave  us  confidence  that  our  covariance  measurements 
were  adequate  measurements  of  the  total  flux  above  the  diffusive  and  wave  boundary  layers  (see 


Veron  et  al.,  2008b). 
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The  instruments  were  deployed  at  the  end  of  the  port  boom  of  RJP  FLIP  approximately  18  m 
from  the  hull  at  an  elevation  of  13  m  above  mean  sea  level  (MSL). 

The  infrared  optical  system  was  set-up  with  the  view  ports  clear  of  the  end  of  the  boom.  The 
meteorological  package  was  placed  15.5  m  from  the  hull  on  the  port  boom  with  all  instruments 
facing  upwind  with  the  exception  of  the  net  radiometer  which  was  deployed  with  its  axis  down¬ 
wind  (figure  1). 

We  also  deployed  2  fast  response,  high  resolution,  subsurface  thermistors  (RBR  ltd  -1040 
95  ms  -  1  Hz  sampling  rate)  placed  at  1.2  m  and  2  m  from  the  mean  sea  level  and  fixed  to 
the  hull  of  RJP  FLIP.  An  upward  looking  waves-enabled  ADCP  (RDI  Workhorse  600kHz),  was 
also  rigidly  mounted  to  FLIP  at  15m  depth  and  yielded  directional  wave  spectra  and  significant 
wave  height  for  the  duration  of  the  experiment.  Finally,  GPS  position  and  RJP  FLIP  heading  were 
sampled  at  50  Hz  and  used  to  correct  the  ADCP  data  for  FLIP  motion  and  to  align  all  other 
directions  to  true  north. 

The  experiment  described  here  took  place  from  20-26  August  2003.  During  this  time,  a  13.5  s 
swell  was  propagating  towards  87°  ±  17°.  For  the  duration  of  the  experiment  the  wind  direction 
was  relatively  constant  from  300°  ±  12°.  The  wind  speed  however  picked  up  from  approximately 
1  m  s_1  at  the  beginning  of  the  week  to  10  m  s-1  a  few  days  later  as  a  storm  passed  by  (see 
figure  2a).  Accordingly,  a  wind  sea  developed  with  waves  of  periods  ranging  from  1.7  s  to  6.3  s, 
propagating  towards  124°  ±12°.  Finally,  the  surface  currents  were  propagating  towards  116°  ±21° 
until  25  August  2003  at  which  point  the  current  turned  significantly.  The  keel  of  FLIP  was  facing 
134°  ±  4°  providing  ideal  conditions  with  wind,  swell,  surface  currents,  and  wind-waves  roughly 
aligned  and  with  the  study  area,  at  and  below  the  port  boom,  clear  of  RJP  FLIP’S  turbulent  wake 


in  both  the  air  and  the  water. 
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3.  Surface  waves  and  velocity  measurements 

To  examine  the  coupling  between  the  surface  temperature,  the  surface  turbulent  velocity  fields 
and  the  surface  waves,  the  surface  wave  elevation  and  direction  of  propagation  must  be  measured. 
This  was  achieved  here  by  single  point  laser  altimetry  measurements  of  the  surface  displacement 
combined  with  the  direct  measurements  of  the  two-dimensional  surface  velocity  over  an  area 
encompassing  the  single  point  altimetry  data.  Details  are  presented  below. 

a.  Surface  altimetry 

The  optical  infrared  system  was  equipped  with  a  12  kHz  single  point  laser  altimeter  with  a 
footprint  on  the  surface  of  about  5  cm  diameter  at  this  operating  height.  It  yielded  instantaneous 
surface  height  rj,  at  50  Hz  after  averaging  the  raw  12  kHz  data.  The  altimeter  relies  on  the 
travel  time  measurement  of  a  light  pulse  which  reflects  from  the  water  surface.  The  altimeter 
performs  significantly  better  when  small  scale  surface  roughness  is  present  at  higher  wind  speeds. 
For  example,  at  0.8  m  s-1  (8.4  m  s-1),  10  %  (67  %)  of  the  raw  returns  led  to  a  valid  height 
measurement.  When  necessary,  dropouts  were  corrected  by  interpolating  the  height  measurements 
using  piecewise  cubic  Hermite  polynomials.  The  footprint  of  the  altimeter  was  situated  22  cm 
from  the  center  of  the  infrared  image.  The  altimeter  data  were  used  to  measure  wave  parameters 
(figure  2)  and  the  instantaneous  distance  from  the  imaging  system  to  the  ocean  surface  in  order  to 
calculate  an  instantaneous  image  resolution  for  the  evaluation  of  the  surface  kinematics. 

Figure  2  shows  the  wind  speed  and  direction  along  with  the  spectragram  of  the  surface  ele¬ 
vation,  and  the  significant  wave  height,  for  the  duration  of  the  experiment.  The  surface  elevation 


Fig.  2 


spectragram  (Figure  2b)  shows  the  swell  peak  at  15  —  20  s  period  and  the  wind-wave  peak  which 
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appears  as  the  wind  rises  from  2  to  8  m  s_1,  and  exhibits  a  classical  downshift  in  frequency  as 
the  wind-wave  field  grows.  Figure  2c  also  shows  good  agreement  between  the  significant  wave 
height  measured  with  the  altimeter  and  that  measured  with  the  waves  ADCP. 

b.  Surface  Velocity 

Fig.  3 


We  have  shown  that  it  is  possible  to  calculate  the  two-dimensional  surface  velocity  u,  and  other 
kinematic  variables  from  infrared  surface  temperature  images  (Veron  et  al.,  2008b).  The  velocity 
at  the  surface  can  further  be  projected  onto  the  directions  parallel  and  perpendicular  to  the  wind  or 
wave  directions.  Here,  we  employ  an  algorithm  typically  employed  for  Particle  Image  Velocime- 
try  (PIV)  and  calculate  the  surface  velocity  by  performing  a  running  normalized  cross-correlation 
on  sub-windows  of  surface  temperature  images  separated  by  33  ms.  Each  cross  correlation  yields 
a  local  average  displacement  (over  the  sub- window)  which  is  then  used  in  conjunction  with  the 
instantaneous  image  resolution  and  time  interval  between  the  successive  images  to  estimate  an 
average  surface  horizontal  velocity.  In  most  cases,  we  used  8x8  pixel  sub-images  with  50%  linear 
overlap  and  image  pairs  separated  by  166  ms  yielding  6  velocity  maps  per  seconds,  each  having 
with  a  spatial  resolution  of  approximately  3  cm.  Details  on  the  performance  of  this  technique  can 
be  found  in  Veron  et  al.  (2008b).  The  surface  velocity  field  u  obtained  is  the  sum  of  velocities 
associated  with  the  mean  currents  uc  ,  the  orbital  velocity  of  the  surface  waves  uw  and  the  turbu¬ 
lence  u'.  Figure  3  shows  examples  of  the  turbulent  surface  velocity  field  if  (i  =  1,2),  obtained 
from  the  deviation  from  the  image  mean  velocity,  and  calculated  from  an  infrared  surface  tem¬ 
perature  image  taken  under  various  wind  speeds.  The  velocity  fields  are  shown  with  the  arrows 
and  are  overlaid  on  the  first  temperature  image  from  which  they  were  calculated.  We  show  here 
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the  turbulent  velocity  rather  than  the  total  velocity  for  easier  visual  comparison  with  the  under¬ 
lying  temperature  structures.  The  velocity  field  shows  a  variety  of  phenomena  with  vortices  and 
regions  of  surface  divergence  and  convergence  that  appear  to  be  correlated  with  the  surface  tem¬ 
perature  structure.  For  example,  the  region  of  warmer  surface  temperature  in  figure  3a  appears  to 
be  associated  with  enhanced  surface  divergence. 

The  benefit  of  this  PIV  technique  is  the  ability  to  simultaneously  study  processes  involved  in 
the  evolution  of  the  surface  temperature  and  surface  velocity  fields  including  mean  currents,  wave 
orbital  motions,  and  turbulence.  In  particular,  direct  surface  velocity  measurements,  along  with 
single  point  altimetry,  allow  us  to  measure  surface  wave  directional  spectra,  and  the  spectra  of  the 
surface  kinematics,  as  illustrated  below. 

c.  Surface  wave -directional  spectra 

For  the  scope  of  the  work  presented  here,  the  single  point  altimetry  from  the  laser  altimeter 
and  the  surface  velocity  from  the  passive  temperature  images  were  combined  to  obtain  the  di¬ 
rectional  spectra  of  the  surface  wave  field.  Velocities  at  each  corner  of  the  infrared  image  image 
were  obtained  by  averaging  all  PIV  estimates  contained  in  local  30  cm  x  30  cm  area.  These  four- 
corner  velocities  were  then  combined  with  the  one  point  surface  elevation  measurements  from  the 
laser  altimeter.  Directional  spectra  were  then  computed  using  the  WAFO  toolbox  (Wave  Analy¬ 
sis  for  Fatigue  and  Oceanography,  Lund  University).  The  directional  wave  frequency  spectrum 
can  be  associated  to  the  cross  spectrum  between  any  two  of  the  considered  measured  quantities, 
i.e.  velocities  and  elevation  (details  can  be  found  in  Young  1994).  To  obtain  estimates  of  the 
full  directional  spectrum,  one  must  choose  analysis  methods  capable  of  inverting  the  relation  be- 
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tween  those  quantities.  We  considered  two  techniques,  the  Iterative  Maximum  Likelihood  Method 
(IMLM,  Pawka  1993),  that  corrects  the  directional  spreading  caused  by  the  conventional  MLM 
technique  (Isobe  et  al.,  1984,  Capon,  1969),  and  the  Maximum  Entropy  Method  (MEM,  Lygre  and 
Krogstad  1986).  Those  two  techniques  showed  some  differences  in  the  computed  spectra,  in  par¬ 
ticular  in  the  spreading  of  the  energy  in  the  vicinity  of  the  wave  frequency  peaks,  but  no  significant 
direction  differences  at  the  wave  energy  peaks  for  swell  and  wind  wave  components.  Wave  peak 
directions  and  frequencies  also  show  a  good  match  with  the  ones  obtained  from  a  Hull-mounted 
RDI  Wave  ADCP.  Only  the  peak  frequencies  and  directions  were  used  in  the  subsequent  analysis. 

4.  Surface  kinematics 

In  addition  to  measuring  the  velocity  field  at  0(1)  cm  spatial  resolution  we  can  also  infer 
the  velocity  gradient  tensor  (Rafael  et  al.,  1998).  This  gives  the  velocity  gradient  tensor 
with  the  same  resolution  and  coverage  as  the  velocity  field.  From  there,  one  can  infer  all  the 
two-dimensional  first  derivative  kinematic  fields  including  divergence  D,  vorticity  u,  shear  S,  and 
normal  deformation  N: 


27T  =  D 

OXn 


^10^ 


V°  lJ 
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(  1  0  ^ 
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+  s 


^01^ 


V  1  0  / 


^  0  -1  ^ 
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or  alternatively, 
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du\  dii2 
dxi  0x2 
du2  dui 
dxi  8x2  ’ 
dui  du2 
dx2  dx±  ’ 
dui  du2 
dx\  dx2  ’ 


(2) 

(3) 

(4) 

(5) 


Here,  we  examine  the  spatial  distribution  of  the  kinematics  fields.  Figure  4  shows  the  two- 
dimensional  wavenumber  spectra  for  the  divergence  and  vorticity  which  were  obtained  from  the 
PIV  processing  of  the  temperature  images.  Figure  4  also  shows  the  two-dimensional  wavenum¬ 
ber  spectrum  for  the  surface  temperature.  The  spectra  were  obtained  by  averaging  120  individual 
spectra  of  a  single  velocity  or  temperature  image,  taken  randomly  during  one  minute  of  data.  The 
data  shown  here  were  collected  under  wind  speeds  of  8  m  s'1  at  the  beginning  of  the  wind  event 
at  05:03  UTC,  on  23  Aug  2003,  during  the  development  of  the  wind  sea.  First,  we  note  that 
the  two  dimensional  spectra  for  the  surface  divergence,  vorticity,  and  temperature,  exhibit  some 
anisotropy  indicating  that  there  is  some  degree  of  directionality  in  the  structures  of  these  spatial 
fields.  Secondly,  we  note  that  the  spectrum  of  the  vorticity  (4b)  is  approximately  aligned  with 
that  of  the  temperature  (4c).  This  is  a  consequence  of  the  coupling  between  surface  temperature 
and  velocity.  We  know  from  previous  work  (Veron  and  Melville,  2001)  that  small  scale  Langmuir 
circulations  and  coherent  turbulence  (streaks)  are  present  at  these  scales  over  this  range  of  wind 
speeds.  Even  if  instantaneous  realizations  of  the  surface  temperature  and  velocity  do  not  clearly 
exhibit  such  structures,  time  averages  can  reveal  their  presence.  Here,  the  two-dimensional  spec¬ 
tra  are  averaged  over  one  minute  of  data.  Colder  streaks  in  the  temperature  images  correspond 
to  water  parcels  that  have  been  organized  by  the  underlying  Langmuir  cells  and  have  collected  in 


Fig.  4 
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windrows  at  the  surface.  These  water  parcels  have  been  in  contact  with  the  air  the  longest  and 
therefore  have  had  time  to  exchange  significant  heat,  leading  to  colder  surface  temperatures.  In 
addition,  these  parcels  of  water  also  have  been  exposed  to  the  wind  stress  the  longest  and  there¬ 
fore,  the  “cold”  surface  temperature  streaks  also  correspond  to  surface  “jets”  in  which  the  surface 
velocity  is  larger  than  that  of  the  surrounding  fluid.  This  is  illustrated  in  Leibovich  (1983).  As  a 
consequence,  there  are  two  horizontal  shear  layers  on  each  side  of  the  temperature  streaks  where 
vertical  vorticity  is  enhanced.  Thus,  when  Langmuir  circulations  are  present,  the  surface  exhibits 
streaks  of  enhanced  velocity,  or  surface  jets,  along  with  streaks  of  colder  temperature  and  streaks 
of  vorticity,  all  in  the  same  direction  and  aligned  with  the  wind.  This  is  precisely  what  is  shown  in 
figure  4.  The  reader  is  reminded  that  the  two-dimensional  spectrum  will  show  the  largest  signal  in 
the  direction  of  the  largest  variance,  i.e.  the  direction  perpendicular  to  the  streaks.  For  example, 
in  the  data  shown  in  figure  4,  the  vorticity  (temperature)  spectrum  peaks  at  52°  (43°)  indicating 
that  the  streaks  in  the  vorticity  fields  are  oriented  along  the  direction  322°  —  142°  (313°  —  133°) , 
roughly  aligned  with  the  wind. 

Finally,  we  note  that  the  spectrum  of  the  divergence  is  perpendicular  to  that  of  the  vorticity. 
This  is  also  consistent  with  the  presence  of  Langmuir  circulations.  Indeed,  Langmuir  circula¬ 
tions  are  generated  from  the  positive  feedback  between  the  surface  jets,  vertical  vorticity,  and 
the  shearing  and  stretching  of  the  vertical  vortex  lines  by  the  Stokes  drift  (Craik  1977,  Leibovich 
1977).  This  means  that  the  surface  waves  are  propagating  in  the  downwind  direction,  with  their 
wavenumber  vector  approximately  aligned  in  the  direction  of  the  surface  streaks.  Accordingly, 
regions  of  divergence  and  convergence  associated  with  the  non-linear  surface  wave  field  and  even 
perhaps  breaking  waves,  will  be  aligned  parallel  to  the  wave  fronts,  perpendicular  to  the  jets  of 
colder  surface  temperature  and  enhanced  vertical  vorticity.  Indeed,  for  the  data  shown  in  figure  4, 


February  2008 


F.  Veron,  W.K.  Melville  and  L.  Lenain 


15 


the  spectral  peak  for  the  divergence  is  approximately  in  the  direction  of  the  wind  and  wind  waves 
indicating  that  the  structure  in  the  surface  divergence  field  is  perpendicular  to  the  predominant 
wind  and  wind-wave  directions,  and  also  perpendicular  to  the  streaks  in  the  surface  temperature 
and  surface  vorticity.  Hence,  when  Langmuir  circulations  (or  Langmuir  turbulence)  are  present, 
the  two-dimensional  wavenumber  spectra  for  the  vorticity  and  divergence  are  perpendicular  to 
one  another  with  the  spectrum  for  the  vorticity  parallel  to  that  of  the  surface  temperature.  The 
presence  of  developing  Langmuir  circulations  during  the  early  stage  of  this  wind  event  was  also 
confirmed  by  visual  inspection  of  the  temperature  images  and  visual  observation  of  the  surface 
during  the  experiment. 

Before  proceeding,  it  is  useful  to  consider  conditions  and  constraints  on  the  surface  velocity 
field  that  lead  to  orthogonality  of  the  two-dimensional  spectra  of  the  divergence  and  vorticity  as  is 
observed  in  figure  4.  Since  the  spectra  shown  in  figure  4  were  averaged  over  one  minute  of  data, 
we  start  with  the  horizontal  components  of  the  mean  surface  velocity  field  defined  by  expanding 
the  velocity  about  an  arbitrary  origin  as  follows: 


Ui  =  u0i  +  ctiX  1  +  (3iX  2  +  7 ix\  +  dixl  +  Xi^  ix2  +  O(xf), 


(6) 


where  the  tq  are  approximated  by  quadratic  functions  of  x,,  with  cubic  and  higher  terms  neglected. 
We  can  then  estimate  the  divergence  and  vorticity  respectively: 


D  =  (ai  +  fa)  +07(271  +  X2)  +  x2(252  +  Xi) 


(7) 


CO 


(a2  -  (3\)  +  27(272  -  xi)  +  x2(X2  ~  26 1) 


(8) 


(9) 


When  the  two-dimensional  spectra  of  the  vorticity  and  divergence  are  orthogonal,  the  structures 
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(streaks  in  the  case  of  Langmuir  circulations)  in  the  divergence  and  vorticity  fields  are  also  or¬ 
thogonal.  Consequently,  the  gradient  on  the  velocity  and  divergence  fields  are  orthogonal  as  well. 
In  other  words, 


which  leads  to  the  following  condition, 


dD  doj 
dxi  dxi 


=  0, 


(10) 


2(7172  -  6162)  -  Xi(7i  +  ^i)  +  *2(72  +  62)  =  0. 


(ID 


Providing  that  Xi  7^  0,  this  leads  to: 

=  -7 i,  (12) 

which  in  turns  leads  to  a  surface  velocity  field  of  the  following  hyperbolic  form 

Ui  =  u0i  +  atxi  +  (3iX 2  +  7 i{x\  -  x\)  +  XiX\x2  +  O(xf),  (13) 


We  deduce  here  that  when  the  surface  velocity  field  is  locally  hyperbolic,  which  leads  to  significant 
stretching,  then  the  two-dimensional  wavenumber  spectra  for  the  vorticity  and  divergence  are  also 
orthogonal.  Hyperbolicity  is  a  sufficient  condition  for  the  spectra  to  be  orthogonal,  but  it  is  not 
necessary.  It  is  unclear  whether  or  not  the  presence  of  Langmuir  circulations  leads  to  hyperbolic 
surface  velocity  fields  and,  to  our  knowledge,  this  question  has  not  been  examined.  However,  it 
is  well  established  that  Langmuir  circulation  will  lead  to  significant  stretching,  and  in  particular 
stretching  of  vortex  lines  which  incidentally  leads  to  the  positive  feedback  and  reinforcement  of 
the  circulations.  These  data  raise  interesting  questions  about  the  nature  of  the  surface  velocity 
fields  during  periods  of  growing  wind  seas  and  developing  Langmuir  circulations. 

In  addition,  the  hyperbolic  surface  velocity  field  in  equation  13  above  immediately  leads  to 
jxxjjj-  —  0-  That  is,  the  two-dimensional  wavenumber  spectra  for  the  shear  and  normal  defor¬ 
mation  are  also  orthogonal.  While  a  hyperbolic  velocity  field  is  not  a  necessary  condition  for 
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this  orthogonality,  it  is  easy  to  demonstrate  that  ffr  =  0  implies  —  0,  and  vice  versa, 

consistent  with  the  orthogonality  of  the  shear  and  normal  deformation  fields  shown  here  (Figure 
5). 

As  the  wind  speed  and  wind  wave  field  stabilize  during  23  and  24  August  2003,  we  observe 
that  the  relationships  between  the  two-dimensional  spectra  of  the  vorticity,  divergence,  shear  and 
normal  deformation  show  some  additional  order.  Not  only  are  the  two-dimensional  spectra  per¬ 
pendicular  in  pairs,  that  is  the  spectrum  of  vorticity  (shear)  is  orthogonal  to  that  of  the  divergence 
(normal  deformation),  but  the  spectrum  for  the  vorticity  (divergence)  is  orthogonal  to  that  of  the 
normal  deformation  (shear).  This  is  clearly  demonstrated  in  Figure  5  which  shows  the  direction 
of  the  wind,  waves  and  currents  (a),  along  with  the  directions  of  the  streaks  (i.e.  perpendicular 
to  the  direction  of  the  peak  spectral  density)  for  the  kinematic  variables  (Figure  5b).  This  further 
apparent  coupling  and  ordering  in  the  surface  kinematics  is  presumably  a  consequence  of  the  sta¬ 
bilization  of  the  surface  wind  wave  field  and  the  wind  drift.  Using  the  simple  formulation  above 
for  the  velocity  fields,  this  additional  orthogonality  yields  a  particular  solution  (though  not  the 
only  one  possible)  with  x%  =  ±(— 1)*27*  which  indeed  leads  to  a  simpler,  more  constrained  form 
of  the  local  surface  velocity  fields 


Ui  =  u0i  +  aiX  1  +  /3iX2  +  7 i{x\  -  x22)  ±  (-l)*27jXi£2  +  O(x^).  (14) 

Fig.  5 

Figure  6  shows  two  examples  of  streamline  patterns  corresponding  to  a  hyperbolic  surface  ve¬ 
locity  field  given  by  equation  14  where  at,  /3i,  and  7 j  where  chosen  randomly.  The  spatial  scale  is 
arbitrary.  As  expected,  the  hyperbolic  surface  velocity  field  shows  some  two-dimensional  stretch¬ 
ing  and  regions  where  the  streamlines  converge  in  a  fashion  that  is  reminiscent  of  the  surface 
velocity  field  that  one  would  expect  in  the  presence  of  Langmuir  circulations.  In  particular,  figure 


18 


JOURNAL  OF  PHYSICAL  OCEANOGRAPHY 


Volume  PKM-000 


5b  shows  that  the  streamlines  of  the  mean  surface  velocity  converge  to  form  surface  streaks. 

Later  still  during  the  week,  on  August  26  2003  and  beyond,  the  geometric  relationship  between 
all  kinematic  fields  is  lost  and  no  obvious  coupling  can  be  found  (figure  5b).  Unfortunately,  the 
surface  currents  turned  significantly  during  these  two  days  (figure  5a).  Consequently,  it  is  difficult 
to  unequivocally  tie  the  apparent  “loss  of  order”  and  the  breakdown  of  the  orthogonality  between 
kinematics  fields  to  either  the  drop  in  wind  speed  and  the  relaxation  in  the  surface  stress,  which 
presumably  destroyed  the  Langmuir  circulations,  or  to  the  change  in  the  direction  of  the  surface 
currents,  which  at  times,  oppose  the  wave  direction  and  Stokes  drift. 


Fig.  6 


Together,  the  data  presented  above  illustrate  that  there  are  clear  couplings  that  occur  on  time 
scales  of  a  minute  and  longer  which  clearly  appear  out  of  long  term  spectral  averages.  The  time 
average  used  here  is  similar  to  that  used  in  the  development  of  the  CLII  mechanism,  and  perhaps 
not  surprisingly,  the  data  shows  that  the  surface  kinematic  fields  are  organized  in  a  manner  con¬ 
sistent  with  the  presence  of  Langmuir  circulations.  These  data  highlight  the  need  for  further  study 
of  the  relationship  between  the  surface  kinematics,  surface  waves,  currents,  and  even  perhaps  the 
surface  temperature  and  gas  concentration  (i.e.  the  surface  molecular  layer). 


5.  Modulation  of  turbulence  by  the  waves 

We  have  seen  in  the  preceding  section  that  the  surface  kinematics  is  coupled  with  the  devel¬ 
oping  wind  seas,  on  time  scales  longer  than  approximately  a  minute.  This  raises  the  question  of 
the  coupling  between  the  surface  turbulence,  waves,  and  currents  on  shorter  time  scales,  and  in 
particular  the  coupling  between  turbulence  and  surface  waves. 


As  described  above,  the  footprint  of  the  laser  altimeter  was  located  within  the  infrared  image 
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and  therefore  within  the  surface  velocity  and  surface  kinematic  fields.  This  allows  us  to  examine 
the  modulation  of  the  surface  velocity  and  kinematic  fields  by  the  waves.  From  the  infrared 
images  and  the  resulting  velocity  images,  along  with  the  wave  directional  spectra  obtained  from 
the  processing  described  in  section  3.c,  we  have  decomposed  the  surface  velocity  field  into  two 
components  u\  and  w2,  which,  respectively,  are  aligned  with,  and  orthogonal  to,  the  direction  of 
propagation  of  the  wind-waves.  To  avoid  sky  reflectance  and  other  effects  in  the  infrared  images, 
only  night-time  data  were  used.  Also,  as  for  the  data  shown  in  figure  3,  we  have  subtracted  image- 
mean  velocity  from  every  velocity  vector  of  that  particular  PIV  estimate.  This  operation  subtracts 
both  currents  uc  and  wave  orbital  velocity  uw,  and  yields  the  surface  turbulent  velocity  field, 
u'.  From  these  we  can  generate  time  series  of  the  statistics  of  the  turbulent  surface  kinematics 
variables  averaged  over  the  footprint  of  the  infrared  imager. 

For  example,  figure  7  shows  time  series  of  the  variance  of  the  turbulent  surface  velocity  in  the 
wave  direction,  uf,  and  the  surface  displacement,  77,  over  a  120s  record  starting  at  05:00  (UTC)  on 
August  24th,  2003.  There  is  a  clear  visual  correlation  between  the  two  variables  with  significant 
bursts  of  turbulence  at  the  surface  that  correlate  with  the  passage  of  surface  waves.  Also,  high 
turbulent  levels  are  located  at,  or  slightly  behind  the  wave  crests.  Spectral  analysis  makes  this  a 
little  clearer. 

Figure  8a  shows  frequency  spectragrams  of  the  surface  displacement  77  for  the  duration  of  the 
FLIP  experiment  in  August  2003.  The  quality  of  the  near-IR  laser  altimeter  data  is  dependent  on 
the  roughness  of  the  sea  surface  with  drop  outs  increasing  for  lower  wind  speeds.  The  data  shown 
here  were  good  for  frequencies  up  to  1  Hz  and  are  low-pass  filtered  at  that  frequency.  The  surface 
displacement  spectra  denoted  by  Sm  is  typical,  showing  both  wind-wave  and  swell  peaks.  The 
wind-wave  peak  appears  with  the  wind  event  (figure  2),  showing  the  classical  downshift  in  wind- 
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Fig.  7 
Fig.  8 


wave  peak  frequency  with  time.  The  spectra  also  show  a  /-4  slope  above  the  wind-wave  peak. 
Figure  8b  shows  that  the  frequency  spectra  for  the  turbulent  velocity  fluctuations  in  the  direction  of 
the  wave  propagation,  u\,  exhibit  peaks  at  the  local  maxima  of  the  wind-wave  spectra.  Figures  8c 
and  d  show  the  squared  coherence  Cu>  v  and  phase  $u'v,  between  u\  and  rj,  respectively.  It  shows 
a  peak  in  the  squared  coherence  of  approximately  0.2  —  0.3,  rapidly  going  to  zero  at  frequencies 
above  approximately  1  Hz.  The  phase  between  the  turbulence  and  surface  displacement  near 
the  peak  of  the  coherence  is  slightly  negative,  indicating  that  the  maximum  in  turbulence  intensity 
lags  the  maximum  surface  displacement.  This  is  presumably  a  consequence  of  the  enhanced  wind- 
stress  and  turbulence  on  the  windward  side  of  the  wave.  This  is  corroborated  by  identical  results 
for  the  other  kinematic  variables  (not  shown).  Incidentally,  this  is  consistent  with  observations 
of  warmer  skin  temperature  also  on  the  windward  side  of  the  waves  that  results  from  enhanced 
wind  stress  and  the  resulting  turbulence  which  leads  to  the  disruption  of  the  aqueous  viscous  layer 
and  the  thermal  molecular  layer,  i.e  the  cool  skin  (Simpson  and  Paulson,  1979;  Miller  and  Street, 
1978;  Veron  et  al.,  2008a). 

Increased  turbulence  and  vorticity  near  the  wave  peak  could  also  be  a  result  of  the  shearing 
effects  on  the  turbulence  and  the  local  stretching  of  vortex  lines  by  the  passing  surface  waves 
(Teixeira  and  Belcher,  2002).  We  examine  next  the  effect  of  the  waves  on  the  turbulence  over 
short,  wave  resolved,  time  scales. 

We  wish  to  extract  the  wave-coherent  turbulent  quantities.  Following  the  standard  Reynolds 
decomposition  u  =  u  +  u'  where  the  overbar  represents  spatial  averages  over  the  PIV  estimates 
(which  are  distributed  over  the  footprint  of  the  infrared  imager).  The  primes  indicate  fluctuating 
quantities  taken  as  deviations  from  the  mean  which  can  be  further  decomposed  into  wave  coherent, 
u,  and  turbulent  components,  u" .  Here  we  use  a  technique  similar  to  that  presented  in  Veron  et 
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al.,  (2008a)  to  extract  the  averaged  wave  coherent  quantities  We  use  the  coherence  and  phase  of 
each  variable  with  the  surface  displacement,  rj,  due  to  the  waves.  Denoting  the  spectrum  of  u'\  by 
Su'  u'^  we  find  that  the  spectrum  of  u\  coherent  with  the  waves  (i.e.  the  spectrum  of  u \ )  is  given 

by  S fa ;  fa 1  Uy  Y  Cu'^ri- 

It  follows  that  the  variance  of  the  phase-coherent  velocity  in  the  direction  of  the  wave  propa¬ 
gation,  (fti),  is  given  by: 

u\  =  I  Sy^du,  (15) 

Fig.  9 

Figure  9  shows  the  variance  of  the  wave-coherent,  along  wave,  surface  velocity  turbulence 
as  a  function  of  the  wind  speed  (a)  and  the  surface  wave  slope2  (b).  While  noisy,  the  data  show 
that  there  is  some  fraction  of  the  surface  turbulence  that  is  coherent  with  the  surface  waves.  The 
data  show  that  the  variance  of  the  wave  coherent  surface  turbulence  increases  with  wind  speed 
and  wave  slope,  with  perhaps  a  tighter  fit  with  the  wave  slope.  Figure  9  also  shows  data  taken  in 
shallow  water  and  at  lower  wind  speeds  from  Scripps  Pier  during  the  winter  of  2003  (see  Veron 
et  al.,  2008b  for  details).  The  two  data  sets  merge  continuously  and  the  large  symbols  are  the  bin 
averaged  data.  The  solid  lines  show  the  linear  fit  through  the  data  (forcing  a  zero  intercept  for 
figure  9b). 

Of  further  interest  is  the  total  fraction  of  this  wave  coherent  turbulence  to  the  total  turbulent 
intensities.  Following  the  work  of  Towsend  (1976,  pp71-77)  and  Texiera  and  Belcher  (2002), 
we  use  the  frame  work  of  rapid  distortion  estimate  this  ratio.  Assuming  initially  isotropic  turbu- 

2Since  the  surface  slope,  {rjx,  rjy ) ,  was  not  directly  measured  in  the  experiments,  the  root  mean  square  wave  slope 

1 

is  instead  approximated  from  the  wave  elevation  frequency  spectrum  by  S'  =  ^  Jq  Svvk2df^j  ,  where  fp  is  the 
peak  frequency  and  k  the  wavenumber  given  by  the  dispersion  relationship.  Finally,  we  estimate  ak  =  y/2S.  This  is 
valid  for  the  gravity  wave  range  only. 
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lence  and  plane  straining  by  a  unidirectional,  monochromatic  linear  wave  field  of  amplitude  a  and 
wavenumber  k,  then 


4  03-r1) 


-n2 


(16) 


u?  5  (/3  +  /3-1)  35 

where,  f3  denotes  the  strain  ratio,  i.e.  the  ratio  of  the  velocity  gradient  tensor  (or  vorticity  compo 
nents)  with  and  without  distortion.  In  this  case,  (3  varies  as: 


(3  «  1  —  ak. 


(17) 


Negelcting  terms  higher  than  quadratic  in  the  wave  slope  ak,  we  find  that: 


u 


4 
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5afc +  -(<*)' 


(18) 


which  indicates  that  the  normalized  wave-modulated  turbulence  in  the  direction  of  the  wave  prop¬ 
agation  is  of  order  0(ak).  It  should  be  noted  that  this  theoretical  estimate  relies  on  assumptions 
of  isotropy  and  homogeneity  in  the  turbulence.  These  can  be  violated  near  the  surface  (Belcher  et 
al.,  1994)  and  in  the  case  where  the  turbulence  in  injected  by  localized  breaking  events;  although, 
as  noted  by  Texiera  and  Belcher  (2002),  even  in  the  latter  case,  the  integral  length  scale  of  the 
turbulence  is  small  compared  to  the  wavelength  of  the  surface  breaking  wave,  hence  satisfying  the 
scale  separation  necessary  for  rapid  distortion  theory  to  be  applicable.  On  that  note,  we  do  not 
expect  the  swell  to  significantly  affect  the  surface  turbulence,  at  least  in  a  rapid  distortion  sense, 
since  the  strain  rate  from  the  orbital  motion  of  the  swell  is  generally  too  weak  compared  to  that  of 
the  turbulence  on  itself.  In  fact,  figure  8c  shows  no  coherence  between  the  waves  and  the  surface 
turbulence  at  the  swell  frequency.  Finally,  the  estimate  above  is  developed  for  monochromatic 


linear  waves  and  not  for  a  full  wave  spectrum,  thus  we  expect  differences  between  this  idealized 
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theoretical  prediction  and  the  observed  field  data.  Nevertheless,  we  feel  that  it  is  valuable  and  in¬ 
teresting  to  compare  the  data  with  the  available  theoretical  estimates.  Figure  10  shows  the  ratio  of 
the  variance  of  the  wave  coherent  turbulence  (uf)  to  the  total  variance  in  the  turbulence  (uf)  as  a 
function  of  the  rms  wave  slope.  The  variance  ratio  in  the  direction  of  the  waves  increases  with  the 
wave  slope.  We  should  note  here  that  we  expect  this  ratio  to  be  relatively  noisy,  especially  at  low 
wind  speeds  and  wave  slopes  as  uf  becomes  smaller.  In  figure  10,  we  also  show  the  theoretical 
estimate  from  equation  18.  While  there  is  a  large  scatter  in  the  data,  and  an  expected  flattening 
due  to  noisy  data  for  low  surface  wave  slopes,  the  theoretical  estimate  agrees  reasonably  well 
with  the  data  at  the  larger  slopes.  These  data  appear  to  be  the  first  field  observations  qualitatively 
supporting  the  use  of  rapid  distortion  theory  to  predict  the  modulation  of  surface  turbulence  by 
surface  waves. 


Fig.  10 


6.  Conclusions 


We  have  shown  that  there  are  periods  of  strong  coupling  between  the  surface  wave  field  and 
wind  forcing,  and  both  the  surface  kinematic  and  temperature  fields.  During  the  early  phase  of  the 
wind  event  described  here,  the  two-dimensional  wavenumber  spectra  of  the  surface  vorticity  and 
surface  temperature,  when  averaged  over  time  scales  of  one  minute  and  longer  show  that  these 
mean  fields  are  organized  in  longitudinal  structures  aligned  with  the  wind  direction.  In  addition, 
we  also  show  that  the  vorticity  and  divergence  structures  at  the  surface  are  orthogonal  to  each 
other.  This  coupling,  which  is  revealed  when  the  surface  fields  are  averaged  over  many  wave 
periods  as  is  done  with  the  CLII  approach,  is  consistent  with  the  presence  of  Langmuir  circula¬ 
tions  (or  Langmuir  turbulence).  As  the  wind  stabilizes,  additional  order  in  the  surface  kinematic 
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fields  is  found  and  the  structures  in  the  normal  deformation  and  shear  are  aligned  with  that  in  the 
fields  of  surface  divergence  and  vorticity,  respectively.  This  orthogonatily  can  be  explained  by  a 
surface  velocity  field  of  hyperbolic  form.  In  turn,  hyperbolicity  leads  to  significant  stretching  and 
convergence  in  the  mean  surface  velocity  field,  once  again  consistent  with  the  effects  of  Langmuir 
circulations. 

Over  short  time  scales  where  the  waves  are  resolved,  we  observe  the  modulation  of  the  surface 
turbulence  by  the  surface  waves.  We  show  that  there  is  a  maximum  wave-coherent  turbulent 
intensity  near  the  crests  of  the  surface  waves.  These  results  are  consistent  with  the  modulation 
of  the  skin  temperature  by  the  surface  waves  where  warmer  temperatures,  resulting  presumably 
from  the  destruction  of  the  molecular  thermal  layer  by  the  turbulence,  occur  near  the  crest  of  the 
waves.  Our  data  shows  that  the  relative  wave-coherent  turbulent  intensities  vary  with  both  wind 
speed  and  surface  wave  slope.  The  data  are  qualitatively  consistent  with  rapid  distortion  theory 
but  more  work  is  needed  for  more  complete  quantitative  comparisons. 
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Figure  Captions 


Fig.  1.  (a)  Instrumentation  setup  from  RJP  FLIP,  (b)  Map  of  the  deployment  area;  RJP  FLIP  was 
moored  off  the  coast  of  southern  California  in  August  2003. 

Fig.  2.  (a)  Wind  speed  (black)  and  direction  (gray)  measured  with  the  sonic  anemometer  for  the 
duration  of  the  experiment,  (b)  Spectragram  of  the  surface  elevation  Sm,  measured  with  the  laser 
altimeter.  Note  the  two  peaks  at  the  frequencies  of  the  swell  and  the  wind  waves,  (c)  Significant 
wave  height  measured  with  the  altimeter  and  compared  with  that  obtained  from  the  RDI  waves 
ADCP  (•). 

Fig.  3.  Infrared  image  of  the  surface  temperature  and  surface  turbulent  velocity  field  calculated 
from  the  passive  temperature  pattern  displacements  between  this  temperature  image  and  one  taken 
33  ms  later,  (a)  Shows  the  temperature  and  surface  velocity  for  a  wind  speed  of  Uio  =  0.5  ms-1 
with  large  up  welling/convection,  and  (b)  for  a  wind  speed  of  Uw  =  3.8  ms-1  where  cold  streaks 
(Langmuir  circulations)  are  clearly  visible. 

Fig.  4.  Two-dimensional  wavenumber  spectra  of  the  surface  divergence  (a),  vorticity  (b)  and 
surface  temperature  (c)  for  a  wind  speed  of  8  m  s-1.  The  spectra  were  calculated  from  an  average 
of  120  individual  spectra  taken  over  one  minute  of  data  starting  at  05:03  UTC  on  23  Aug  2003. 
(d)  Shows  the  directions  for  the  wind,  from  305°,  the  wind  waves,  to  117°,  the  swell,  to  97°,  the 
vorticity  spectral  peak,  52°,  the  divergence  spectral  peak,  122°,  and  the  temperature  spectral  peak, 


43°. 
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Fig.  5.  (a)  Direction  of  the  wind,  wind-waves,  swell  and  current,  (b)  Directions  of  the  streaks  (i.e. 
the  direction  perpendicular  to  that  of  the  peak  spectral  densities  and  the  maximum  variance),  for 
the  vorticity,  divergence,  shear,  and  normal  deformation. 


Fig.  6.  Examples  of  streamline  patterns  showing  hyperbolic  surface  velocity  fields  as  in  equation 
14  where  cu,  and  7*  where  chosen  randomly.  The  scale  is  arbitrary. 


Fig.  7.  Time  series  of  surface  displacement  and  the  spatial  variance  of  the  along  wave  surface 
turbulence  over  a  120  s  record  starting  at  05:00  (UTC)  on  August  24th,  2003. 


Fig.  8.  Spectragrams  of  (a)  the  surface  displacement  7  and  (b)  the  along  wave  turbulent  surface 
velocity  u\,  at  night,  (c)  and  (d)  show  the  squared  coherence  and  phase  between  the  surface 
displacement  and  the  turbulent  surface  velocity,  respectively  . 


Fig.  9.  Variance  of  the  wave-coherent,  along  wave,  surface  velocity  turbulence  with  the  wind 
speed  (a)  and  the  rms  surface  wave  slope  (b).  The  large  symbols  show  the  bin-averaged  data  and 
the  solid  lines  are  the  fits  to  the  bin  averaged  data. 


Fig.  10.  Fraction  of  variance  of  the  wave  coherent  turbulence  (ul)  to  the  total  variance  in  the 
turbulence  (uf  )  as  a  function  of  the  wave  slope.  The  large  symbols  show  the  bin-averaged  data 
and  the  solid  gray  line  shows  the  rapid  distortion  theory  prediction  of  equation  18. 
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Fig.  1.  (a)  Instrumentation  setup  from  R/P  FLIP,  (b)  Map  of  the  deployment  area;  R/P  FLIP  was 
moored  off  the  coast  of  southern  California  in  August  2003. 
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Fig.  2.  (a)  Wind  speed  (black)  and  direction  (gray)  measured  with  the  sonic  anemometer  for  the 
duration  of  the  experiment,  (b)  Spectragram  of  the  surface  elevation  Sm,  measured  with  the  laser 
altimeter.  Note  the  two  peaks  at  the  frequencies  of  the  swell  and  the  wind  waves,  (c)  Significant 
wave  height  measured  with  the  altimeter  and  compared  with  that  obtained  from  the  RDI  waves 


ADCP  (•). 
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Fig.  3.  Infrared  image  of  the  surface  temperature  and  surface  turbulent  velocity  field  calculated 
from  the  passive  temperature  pattern  displacements  between  this  temperature  image  and  one  taken 
33  ms  later,  (a)  Shows  the  temperature  and  surface  velocity  for  a  wind  speed  of  Uiq  =  0.5  ms-1 
with  large  up  welling/convection,  and  (b)  for  a  wind  speed  of  Uw  =  3.8  ms-1  where  cold  streaks 
(Langmuir  circulations)  are  clearly  visible. 
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Fig.  4.  Two-dimensional  wavenumber  spectra  of  the  surface  divergence  (a),  vorticity  (b)  and 
surface  temperature  (c)  for  a  wind  speed  of  8  m  s_1.  The  spectra  were  calculated  from  an  average 
of  120  individual  spectra  taken  over  one  minute  of  data  starting  at  05:03  UTC  on  23  Aug  2003. 
(d)  Shows  the  directions  for  the  wind,  from  305°,  the  wind  waves,  to  117°,  the  swell,  to  97°,  the 
vorticity  spectral  peak,  52°,  the  divergence  spectral  peak,  122°,  and  the  temperature  spectral  peak, 


43°. 
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Fig.  5.  (a)  Direction  of  the  wind,  wind-waves,  swell  and  current,  (b)  Directions  of  the  streaks  (i.e. 
the  direction  perpendicular  to  that  of  the  peak  spectral  densities  and  the  maximum  variance),  for 
the  vorticity,  divergence,  shear,  and  normal  deformation. 


Fig.  6.  Examples  of  streamline  patterns  showing  hyperbolic  surface  velocity  fields  as  in  equation 
14  where  on,  and  7,  where  chosen  randomly.  The  scale  is  arbitrary. 
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Fig.  7.  Time  series  of  surface  displacement  and  the  spatial  variance  of  the  along  wave  surface 
turbulence  over  a  120  s  record  starting  at  05:00  (UTC)  on  August  24th,  2003. 
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Fig.  9.  Variance  of  the  wave-coherent,  along  wave,  surface  velocity  turbulence  with  the  wind 
speed  (a)  and  the  rms  surface  wave  slope  (b).  The  large  symbols  show  the  bin-averaged  data  and 
the  solid  lines  are  the  fits  to  the  bin  averaged  data. 
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Fig.  10.  Fraction  of  variance  of  the  wave  coherent  turbulence  (uj)  to  the  total  variance  in  the 
turbulence  (uf)  as  a  function  of  the  wave  slope.  The  large  symbols  show  the  bin-averaged  data 
and  the  solid  gray  line  shows  the  rapid  distortion  theory  prediction  of  equation  18. 


